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ABSTRACT 

This paper presents a composite diagnosis of synoptic-scale forcing mechanisms associated with extratropical 
cyclone evolution. Drawn from 12 cyclone cases that occurred over the continental United States during the cool 
season months, the diagnosis provides a "climatology" of development mechanisms for difference categories 
of cyclone evolution ranging from cyclone weakening through three stages of cyclone intensification. Compu- 
tational results were obtained using an "extended" form of the Zwack-Okossi equation applied to routine 
upper-air and surface data analyzed on a 230 km x 230 km grid. Results show that cyclonic vorticity advect.on 
which maximizes in the upper troposphere, was the primary contributor to cyclone development regardless of 
the sta°e of development. A second consistent contributor to development was latent heat release. Horizontal 
temperature advection. often acknowledged as a development mechanism, was found to contribute to develop- 
ment only during more intense stages. During weakening and weaker development stages, temperature advect.on 
opposed development, as the warm-air advection invariably found at upper levels was dominated by cold air 
advection in the lower half of the troposphere. In the more intense stages, development was moderated by dry- 
adiabatic cooling associated with the ascending vertical motions. 


1. Introduction 

Studies of extratropical cyclone evolution have 
been for much of this century a major interest of the 
meteorological community, largely because they are 
the source of much of the inclement weather found 
in the midlatitudes and because, with their anticy- 
clone partners, they are responsible for much of the 
potential to kinetic energy transformation found in 
the general circulation (Saltzman 1970). Dating 
from the work of the “Bergen school in the 1910s 
and 1920s (e.g„ Bjerknes 1919; Bjerknes and Sol- 
berg 1922), through the pioneering work of Sut- 
cliffe (e.g., Sutcliffe 1939, 1947), and finally into 
contemporary interest in “explosive” cyclogenesis 
(e.g., Sanders and Gyakum 1980; Roebber 1984; 
Sanders 1986), a great deal has been learned about 
extratropical cyclone evolution. In some cases con- 
firming previous understanding and in other cases 
providing new insight, recent papers demonstrate 
that cyclone development is often associated with 
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cyclonic vorticity advection aloft (Sanders 1986, 
MacDonald and Reiter 1988); upper-tropospheric/ 
lower-stratospheric warm-air advection (Hirsch- 
berg and Fritsch 1991a,b; Lupo et al. 1992); latent 
heat release (Gyakum 1983a, b; Kuo and Reed 
1988); surface energy fluxes (Nuss and Anthes 
1987; Kuo et al. 1991); lowered static stability 
(Nuss and Anthes 1987; Smith and Tsou 1988); 
stratospheric extrusion and subsequent advection 
of high-potential-vorticity air (Uccellini et al. 
1985; Elsberry and Kirchoffer 1988); and fa- 
vorable jet streak positioning (Uccellini et al. 
1984; Uccellini and Kocin 1987). Comprehen- 
sive reviews of cyclone studies can be found in 
Hoskins (1990), Uccellini (1990), and Bosart 
(1994). 

This large array of antecedent studies has pro- 
vided remarkable insight into cyclone development 
processes. However, while the basic cyclogenetic 
mechanisms may be well understood, the “mix” of 
these mechanisms during different stages of cyclone 
development is not always so clear. With this in 
mind, this paper describes a composite study of a 
set of forcing mechanisms that is sufficiently com- 
plete to account for the bulk of the synoptic-scale 
cyclone development process. The composites are 
drawn from 12 extratropical cyclones that occurred 
from late fall to early spring over the continental 
United States and extreme southern Canada. Using 
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pressure tendencies diagnosed at the cyclone center 
to represent cyclone development. 54 individual 
map times represented by these cyclones were strat- 
ified (sample size in parenthesis) into weakening 
(8). weak development (9), moderate development 
(26). and strong development (11) categories. Di- 
agnostic fields were then averaged for each category 
to provide composite values of forcing processes 
representative of different stages of cool-season cy- 
clone evolution. 


2. The Zwack-Okossi equation 

The primary diagnostic tool used in this study is the 
Zwack-Okossi (Z-O) equation, originally derived in 
quasigeostrophic form by Zwack and Okossi (1986) 
and recently redeveloped in a generalized form by 
Lupo et al. (1992). Combining the vorticity equation 
and a form of the first law of thermodynamics through 
the hydrostatic equation, the generalized Z-0 equa- 
tion, given as Eq. (9) in Lupo et al. (1992), is 

k*V X F dp 
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where p, is a near-surface pressure level, p, is an upper 
pressure level sufficiently high to encompass most of 
the atmosphere, is the near-surface geostrophic vor- 
ticity, is the absolute vorticity (C, +/), / is the Corio- 
lis parameter, R is the dry-air gas constant, V is the 
horizontal wind vector with components u and u, £> is 
the diabatic heating /cooling rate per unit mass, c p is the 
specific heat at constant pressure, S is the static stability 
parameter [ ( TI8){d9ldp ), where 6 is potential tem- 
perature], w is the vertical motion in isobaric coordi- 
nates ( u> = dpldt), F is the frictional force, and PD 
= (.Pi ~ P .)~ x • 

Lupo et al. (1992) suggest that when synoptic-scale 
or larger phenomena are the object of the diagnosis it 
is appropriate to use the following simplified form of 
( 1 ) 



(A) 


— V • VT + j +Swjyjdp, (2) 

(B) (C) (D) 

the so-called “extended 1 ’ form of the Z-0 equation. 
This simplified form is obtained from ( 1 ) by eliminat- 
ing terms F, G, and H based on scaling arguments sim- 
ilar to those applied to the height tendency equation of 
Tsou et al. (1987) and term I because it integrates to 
near zero as a result of compensating convergence/di- 
vergence. The frictional term (E) was eliminated be- 
cause it is normally a small contributor and because its 
computation requires making assumptions about sur- 
face conditions that are not known to any degree of 


certainty. Note that nonquasigeostrophic forcing is in- 
cluded in the right-hand-side terms. Also note that (2) 
is in a sufficiently simple form that the calculations for 
a particular map time can be accomplished quickly 
enough for near real-time diagnosis of cyclone events. 
Although not important for this paper, this feature has 
proven useful in the utilization of tine Z-O methodol- 
ogy in the classroom. 

The reader may wonder why the Z-0 equation is 
used rather than other alternatives. Equation (2) is a 
generalization of the Petterssen-Sutcliffe equation 
(Petterssen 1956, 320-325) that permits explicit con- 
sideration of both dynamic and thermal forcing mech- 
anisms at all atmospheric levels and includes nonqua- 
sigeostrophic forcing. Alternative equations that might 
be used for the diagnosis presented here either do not 
explicitly include forcing at all levels (e.g., Petterssen- 
Sutcliffe equation), do not explicitly isolate dynamic 
forcing mechanisms (e.g., Hirschberg-Fritsch equa- 
tion; Hirschberg and Fritsch 1 99 1 a,b) , or do not lend 
themselves to easy isolation of the forcing contribu- 
tions from each level (e.g., height tendency equation; 
Tsou et al. 1987). The integrands of the p t to p t inte- 
gration are proportional to the horizontal divergence 
forced at each pressure level p by horizontal vorticity 
advection (vadv, term A), horizontal temperature ad- 
vection (tadv, term B), diabatic heating (only latent 
heat release in this study, lath, term C), and adiabatic 
cooling/ warming (adia, term D). 

Term A represents the change in near-surface geo- 
strophic vorticity forced by divergence/convergence 
induced by ageostrophic motions that occur as the at- 
mosphere attempts to reestablish geostrophic balance 
in response to local horizontal vorticity advection. 
Thus, for example, cyclonic vorticity advection at any 
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Fig. 1. Total analysis domain represented by all grid points. 
Computational (display) domain enclosed by dashed (solid) 
boundaries. 


level above the surface will induce horizontal diver- 
gence at that level, which in turn will contribute to geo- 
strophic vorticity (pressure) increases (decreases) at 
the surface. Terms B and C represent the horizontal 
divergence/convergence and change in near-surface 
geostrophic vorticity forced by integrated horizontal 
temperature advection and diabatic heating/cooling, 
respectively, from the surface to level p (interior inte- 
gral of these terms). For example, horizontally differ- 
ential (given by the Laplacian operator), integrated 
warming from p , to p will elevate the pressure levels at 
level p and above. This in turn will force ageostrophic 
outflow (divergence) and surface vorticity (pressure) 
increases (decreases). The examples given for terms 
A-C also typically yield upward vertical motions. 
Such motions force dry adiabatic cooling, which in turn 
for a stable atmosphere (5 > 0) yields a value for term 
D that is opposite to that of terms A-C. Thus, term D 
is often found to oppose the development characteris- 
tics of terms A-C. 

3. Computational features 

a. Data, analyses, and finite differencing 

Routinely observed upper-air and surface data over 
the continental United States and surrounding regions 
were used in this study. The upper-air data were ob- 
tained from standard rawinsonde soundings taken twice 
daily at 0000 and 1200 UTC. From the sounding data, 
the surface pressure (mb) and the geopotential heights 
(m), temperatures (C), w'ind speeds (kt), and wind 
directions (°) at the mandatory reporting levels up to 
50 mb were retained for computations. From the sur- 
face data, sea level pressures (mb) at 0000, 0600, 1200, 
and 1 800 UTC and 6-h precipitation accumulations [in. 
(6 h) _I ] recorded at 0000 and 1200 UTC were re- 
tained. All data were obtained using the Purdue Uni- 
versity Weather Processor (WXP) , which has a routine 
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Fig. 2. Response function for Barnes objective analysis. 


to sort and format the data arriving through the do- 
mestic data circuit. 

All surface and mandatory level data were fit to a 25 
X 17 polar-stereographic grid (Fig. 1 ) using a two-pass 
Barnes objective analysis scheme (Barnes 1964, 
1973). The polar stereographic grid is true at 60°N with 
a grid spacing of 230 km X 230 km. Application of the 
Bames objective analysis requires the specification of 
two filter parameters, C and G, which were chosen as 
60 000 and 0.8, respectively. This provides a response 
function (Fig. 2) in which essentially no information 
is retained at wavelengths less than 500 km, approxi- 
mately 23% is retained at 1000-km wavelengths, and 
83% is retained at 2000-km wavelengths. All calcula- 
tions involved with gridding the data were performed 
through the algorithms contained in WXP. After the 
grids at the mandatory levels were created, these data 
were then interpolated linear in ln(p) to standard iso- 
baric surfaces to form a set of data vertically distributed 
in 50-mb increments from 1050 to 50 mb. In perform- 
ing the interpolation, ail winds below ground were set 
to zero, and geopotential heights were calculated with 
the hypsometric equation using the interpolated tem- 
peratures. All horizontal derivative calculations at grid 
points along and within the interior solid box in Fig. 1 
were determined with fourth-order finite differencing, 


SHAPIRO FLTER RESPONSE FUNCTION 



Fig. 3. Response function for Shapiro filtering procedure. 
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while those along the dashed lines were determined 
with second-order finite differencing. Second-order fi- 
nite differencing was also applied to all vertical deriv- 
atives. Vertical integrals were calculated using the trap- 
ezoidal rule. 


b. Vertical velocities 


Vertical velocities were calculated using the omega 
equation in the “extended” form as given by Tsou et 
al. (1987), 


f 2 d 3 

v! + 7« +/, «?j 


" = -r< v ' v W 

<7 Op 


-f — V 2 (V-V7)- — V 2 (-V (3) 
T op op \c p J 


where o is the static stability parameter [a = -(/?77 
P9)(d9/dp)]. Sequential overrelaxation was used to 
solve (3) with the relaxation iteration being terminated 
when the absolute difference between two successive 
solutions of uj was less than 10 " 3 pb “ 1 at all grid points. 
In the relaxation, the initial condition, the upper vertical 
boundary, and the lateral boundary conditions were set 
to zero. The lower vertical boundary value was speci- 
fied as —V * Vpsgz i where z is the height of the terrain 
and p s is the surface air density calculated from the 
surface pressure and temperature. Vertical velocities 
for grid points below the earth’s surface were assigned 
a value of zero. 


c. Diabatic heating 

In this study, the only diabatic heating contribution 
included was latent heat release (LHR). The LHR was 
calculated by taking the surface 6-h precipitation ac- 
cumulation ending at each computational map time 
(0000 or 1200 UTC) and converting it to its equivalent 
integrated latent heating rate. The relationship relating 
precipitation rate to latent heating rate can be derived 
through the basic relationship 

<4> 

where £) is the latent heating rate per unit mass, L is 
the latent heat of vaporization, and dwjdt is the rate 
of change of the saturation mixing ratio. From (4), the 
total column heating rate (H) can be calculated by 

(5) 

where A p is the depth of the column in which the LHR 
is occurring. The precipitation rate P can be expressed 
as 

<6) 


where p is the air density. Applying the hydrostatic 
equation to (6) gives 



(7) 


where g is the acceleration due to gravity (9.81 m s ! ). 
Finally, substitution of (7) into (5) yields 



( 8 ) 


The relationship given in (8) provides the means to 
convert a precipitation rate to an equivalent latent heat- 
ing rate. This heating rate was then distributed in the 
vertical in the following manner: 

(i) if the gridpoint surface pressure was 1000 mb or 
greater, LHR was bounded (set equal to zero) at the 
lower and upper levels by 900 and 300 mb, respec- 
tively, with the maximum LHR at 700 mb; 

(ii) if the gridpoint surface pressure was less than 
1000 mb but greater than or equal to 900 mb, LHR was 
bounded by 850 and 300 mb with the maximum LHR 
at 700 mb; 

(iii) if the gridpoint surface pressure was less than 
900 mb but greater than or equal to 850 mb, LHR was 
bounded by 800 and 300 mb with the maximum LHR 
at 650 mb; and 

(iv) if the gridpoint surface pressure was less than 
850 mb, LHR was bounded by 750 and 300 mb with 
the maximum LHR at 600 mb. 


With this distribution method, no matter what the 
surface elevation, LHR began at approximately 50- 
150 mb above the surface with the maximum LHR 
placed approximately 200-300 mb above the sur- 
face. The location of the LHR maximum is consis- 
tent with the mean latent heat release profiles re- 
ported for different cyclone intensities and devel- 
opment phases by MacDonald and Reiter (1988). 
Once the appropriate distribution method from 
above was chosen for a grid point, the LHR was then 
increased linearly in In (p) from the lower bound to 
the level of maximum LHR and then decreased lin- 
early in ln(p) from the level of maximum LHR to 
the upper bound. 

To examine the sensitivity of the results to the 
vertical distribution of LHR, two alternative vertical 
distributions were tried for one map time, 1200 UTC 
3 December 1990, a time of significant LHR influ- 
ence. One alternative was to place the maximum 
LHR 50 mb above the lower bound (approximately 
100-200 mb above the surface), while the other 
was to place the maximum 400 mb above the lower 
bound (approximately 450 to 550 mb above the sur- 
face). The upper and lower bounds were the same 
as determined in the original procedure. Consider- 
ing the impact on both terms C and D, the pressure 
tendency at the cyclone center was 8.5% larger for 
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T*C (utC) ( a ) 


FlG. 4. (a) Central sea level pressures (mb) and (b) cyclone tracks for the 12 cyclone cases. 


the lower maximum and 23% smaller for the higher 
maximum. Further, maps of pressure tendency (not 
shown) showed little change in the pattern and only 
small differences in magnitudes at most grid points. 
Thus, we conclude that the basic conclusions of this 
paper are not significantly influenced by the choice 
of LHR maximum level. 


d. Pressure tendency calculations 

The vertical computational domain extended from a 
near-surface level ( pt ), defined as the first pressure 
level above the surface identically divisible by 50 (950 
or 900 mb at most grid points) to 50 mb (p,). Forcing 
quantities in (2) were calculated at each level and then 
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Fig. 4. ( Continued ) 


integrated to yield geostrophic vorticity tendencies at 
level p t . At this level, the tendencies are always above 
the surface but are considered close enough (from 0 to 
49 mb above the surface) to be representative of actual 
surface development. After the geostrophic vorticity 
tendency field was calculated for each term in (2), each 
term was filtered using a two-dimensional, fourth-order 


Shapiro (1970) filtering scheme. Filtering was neces- 
sary to remove subsynoptic “noise” introduced 
through finite differencing in the basic calculations. 
The response function for this filter (Fig. 3) retains 
essentially zero information at wavelengths less than 
1000 km, approximately 40% at 1500 km, and 90% at 
2000 km. Prior to filtering, the gridpoint values along 
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Cyclon# Track 




the dashed lines shown in Fig. 1 were used to expand 
the grid by two rows at the top and bottom and two 
columns on each side. For example, the values at the 
grid points along the top row dashed line were used as 
die values for the grid points in the new row above the 
dashed line. This procedure was performed on all sides 


of the grid using the gridpoint values along the other 
dashed lines. Each of the four comer points was set to 
the gridpoint value diagonal from it. After this first ex- 
pansion, the procedure was performed once more using 
the values at the grid points determined from the first 
expansion. Expanding the grid was necessary in order 
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Fig. 5. Distribution of correlation coefficients between 2-0 
and observed pressure tendencies for all map times. 


Fig. 6. The Z-0 versus observed pressure tendencies at the grid 
point of the cyclone center for all map times. The equation for simple 
linear regression is provided on the plot. 


to reduce the boundary effects, which can pose prob- 
lems with small grids like the one used in this study. 

The filtered geostrophic vorticity tendencies were 
summed to provide the total Z-O geostrophic vorticity 
tendency field at p t . Then, sequential overrelaxation 
was used to obtain the total Z-0 near-surface height 
tendency field. In this relaxation, the initial condition 
was set to zero, while the lateral boundaries were set 
to the observed 6-h height tendencies, computed from 
the present minus previous 6-h sea level pressure val- 
ues using Az/A t = ( 1 /p^)( A/? SL / Af), where p s l is 
the sea level pressure. Relaxation continued until the 
absolute difference between two successive solutions 
of height tendency was less than 10' 5 m s _1 at all grid 
points. The final Z-O height tendency field was then 
converted to a surface pressure p s tendency field using 
{dpjdt) = p s g(dztfd f). 

e. Compositing procedure 

Compositing was done by averaging Z-0 near-sur- 
face geostrophic vorticity and surface pressure tenden- 
cies and vertical profiles of relevant parameters at the 
cyclone center (the grid point with lowest sea level 
pressure) for each development category. The center 
point was used because it is an easily identified point 
common to all cyclones and because at this point the 
propagation component of the pressure tendency is 
equal to zero. Thus, the center point pressure tendency 
corresponds to the quasi-Lagrangian pressure change 
or to the development rate. 

Individual map times of cyclone development were 
stratified into four development categories based on the 
center-point Z-0 surface pressure tendencies at each 
time. The four categories are defined as weakening 
(WK, dpjdt > 0); weak development [WD, 0 ^ dpj 
dt> -4 mb (12 h)" 1 ]; moderate development [MD, 
-4 mb (12 h)~ l ^ dpjdt > -8 mb (12 h)" 1 ] (ap- 
proximately 0.5-0.9 bergerons); and strong develop- 


ment [SD, -8 mb ^ dpjdt (12 h) _1 ] (approximately 
0.9 bergerons and greater). Using this classification 
method, 8, 9, 26, and 1 1 map times were included in 
the WK, WD, MD, and SD groups, respectively. Re- 
sults are presented as center-point averages over each 
development category. 

4. Cyclone cases 

Examined in this study were the 12 cyclones iden- 
tified by the cyclone tracks and central pressures shown 
in Fig. 4. These cases were chosen to provide a repre- 
sentative sample of typical cool season cyclones that 
included examples of all of the development categories 
and were located in the interior of the data network. 
Examination of the central sea level pressure of each 
cyclone indicates phases of development ranging from 
weakening to explosive intensification. Four cy- 
clones — cases 03DEC90, 02NOV92, 12NOV93, and 
24MAR94 — clearly exhibited both a period of deep- 
ening and a period of weakening. Two cyclones — 
cases 12NOV93 and 26JAN94 — experienced brief pe- 
riods (6 h) of explosive development at rates of 1.3 
and 1.1 bergerons, respectively, while case 24MAR94 
developed explosively for 18 h at a rate of 1.1 berge- 


Table 1. Mean absolute values [MAV, mb (12 h)' 1 ] and 
correlation coefficients (COR) for the four development classes. 
DIFF represents the percentage difference between OBS and Z-0 
MAV values, that is, [(Z-0 MAV) — (OBS MAV)]/(OBS MAV) 
X 100. 



WK 

WD 

MD 

SD 

OBS MAV 

4.20 

5.24 

4.51 

6.02 

Z-0 MAV 

3.55 

3.84 

3.99 

6.05 

DIFF 

-15 

-27 

-12 

0,5 

COR 

0.77 

0.81 

0.82 

0.83 
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FiO. 7. Composite near-surface geostrophic vorticity tendencies <10~* *-) at the cyclone center forced by the vort city 
advection term (vadv). temperature advection term (tadv). adiabatic term (adia). and latent heating term (lath) he total 
Z-0 near-surface geostrophic vorticity tendency (zoto); and the Z-O surface pressure tendency [zopt. mb (12 h) J 
for each of the four development categories. 


rons. Only one case, 12NOV92, experienced an ex- 
tended explosive development period that lasted at least 
24 h. During the period 0600 UTC 12 November-0600 
UTC 13 November, the cyclone intensified at a rate of 
1.9 bergerons. The other cyclone cases exhibited vari- 
ous degrees of nonexplosive deepening with some pe- 
riods of weakening. 

5. Results 

a. Comparisons of Z-O and observed pressure 
tendencies 

The Z-O surface pressure tendencies were com- 
pared with observed (OBS) sea level pressure changes, 
computed from observed sea level pressure at the di- 
agnosis time minus observed sea level pressure 6 h ear- 
lier. Adoption of the backward 6-h time difference was 
dictated by classroom use of this methodology for real- 
time cases. This provides us with an observed tendency 
field that is available for comparison with the Z-O ten- 


dencies as soon as they are calculated. The correlation 
coefficients (COR) for all map times in this study (Fig. 
5) suggest very favorable pattern comparability be- 
tween the two fields. In fact, 34 of the 54 map times 
had COR values exceeding 0.8, and 14 of those values 
were above 0.9. The worst and best COR values were 
0.50 and 0.96, respectively, and the average was 0.81. 
Furthermore, in comparisons of Z-O pressure tenden- 
cies to OBS pressure changes at the gridpoint locations 
of the cyclone center (Fig. 6), there were only five 
occurrences of Z-O and OBS tendencies having op- 
posite sign, and in each of those five, one of the ten- 
dencies was always near zero. This suggests that the 
Z-O equation provided a good diagnosis of the forcing 
processes leading to surface pressure changes, partic- 
ularly near the cyclone center. Closer examination of 
Fig. 6 reveals that the Z-O tendencies were biased to- 
ward values smaller than observed. Finally, the mean 
absolute values (MAV), averaged over all grid points 
and map times, also suggest good comparability be- 
tween the two fields. The MAV for the Z— O pressure 
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Fig. 8. Composite vertical profiles of horizontal absolute vorticity advection (1CT 9 s 2 ) 
at the cyclone center for each of the four development categories. 


tendencies, at 4.31 mb ( 12 h)"\ was only 12% smaller 
than the OBS MAV, at 4.89 mb (12 h)~ l . Thus, while 
the COR values suggest good comparability of pressure 
fall and pressure rise regions between the Z-0 and 
OBS fields, the similar MAV’s suggest good compar- 
ability between the magnitudes in the pressure fall and 
pressure rise regions. 

Table 1 shows the comparison statistics, computed 
as averages, between the Z-0 and OBS pressure ten- 
dencies in each of the four development categories. The 
WK phase had the lowest COR value, while the SD 
phase had the best COR value, as well as the smallest 
difference between OBS and Z-0 MAV values at 
0.5%. The \VD phase had the largest difference in 
MAV values, with the Z-0 being 27% smaller than 
the OBS. In general, the comparisons between the Z- 
O and OBS pressure tendency fields became more fa- 
vorable with stronger development. This suggests that 
the Z-0 model works best when the synoptic-scale 
flow regime is active with prominent, intensifying cy- 
clones. 


b. Composite Z—O tendencies 

Figure 7 shows the composite near-surface geo- 
strophic vorticity tendencies at the cyclone center 
point forced by each term in (2), as well as the total 
near-surface vorticity tendency (zoto) and the Z-0 
surface pressure tendency (zopt) for each develop- 
ment category. As expected, zoto (zopt) was nega- 
tive (positive) for weakening cyclones and increas- 
ingly positive (negative) for increasingly stronger 
cyclone development. 

The most consistent contributor to development was 
the vorticity advection, which monotonically increased 
in magnitude for stronger development categories. Fail- 
ure of cyclones to intensify during the weakening pe- 
riods could be largely attributed to the temperature ad- 
vection influence, which made the largest negative con- 
tribution. This negative contribution also occurred 
during weak development periods and would have re- 
sulted in cyclone weakening except for the added pos- 
itive influence of latent heat release. For both of these 
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Fig. 9. As in Fig. 8 except for horizontal temperature advection (10"* K s' 1 ). 


phases, dry-adiabatic processes had little impact on de- 
velopment. 

During moderate development, the temperature ad- 
vection influence was nearly zero, while the latent heat 
release influence remained nearly the same as during 
weak development. Adiabatic cooling assumed a 
stronger role, but its negative influence was more than 
compensated by the vorticity advection. Finally, strong 
development is marked by increased positive contri- 
butions by the dominant vorticity advection as well as 
temperature advection and latent heat release, which 
together were able to more than offset the increased 
adiabatic cooling influence. This latter result is quali- 
tatively similar to the results of Lupo et al. (1992) for 
cases of marine and land explosive cyclogenesis. 

c. Vertical profiles 

Insight into the physical mechanisms responsible for 
the preceding results is best achieved by examining 
vertical profiles of the basic physical process repre- 


sented in each term, that is, horizontal vorticity advec- 
tion, horizontal temperature advection, latent heat re- 
lease, and adiabatic temperature change (Scu). Com- 
posite profiles for each development category are 
displayed in Figs. 8-13. 

The composite profiles of horizontal absolute vortic- 
ity advection are shown in Fig. 8. The profiles for the 
WD, MD, and SD categories show a cyclonic vorticity 
advection (CVA) maximum in the upper troposphere, 
where the strongest winds are found, with the maxi- 
mum becoming more prominent and well defined with 
increasing intensity of surface development. At low 
levels, weak anticyclonic vorticity advection (AVA) is 
seen for these three composites, with magnitude and 
vertical depth of influence that decreased as surface 
development increased. These profiles are similar to 
those reported in MacDonald and Reiter s ( 1988) com- 
posite study. Also, the SD profile is similar to the vor- 
ticity advection profiles in Lupo et al. ( 1992) . The Wk 
composite showed w'eak vorticity advections at all lev- 
els with no prominent CVA or AVA maximum. How- 
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ever, similar to the other categories, CVA was found 
at all but the lower portion of the troposphere. Note the 
good correspondence between increasing CVA in the 
four categories and the increasing near-surface geo- 
strophic vorticity tendency values contributed by vadv 
in Fig. 7. 

Similar to vorticity advection, temperature advection 
profiles show a warm-air advection ( WAA) maximum 
at upper atmospheric levels (Fig. 9), a common feature 
of cyclone development emphasized by Hirschberg and 
Fritsch (1991b) and confirmed by Lupo et al. (1992). 
Notice that as the degree of surface development in- 
creased in the four composites, the value of the upper- 
level WAA maximum also increased. Also, in the MD 
and SD composites, a secondary WAA maximum was 
located in the lowest 150 to 200 mb, similar to the 
results of Lupo et al. ( 1992). Although WAA increased 
with increased surface development, it became domi- 
nant only in the SD composite, when as shown in Fig. 
7 it forced near-surface geostrophic vorticity increases. 
In all four classes, cold-air advection (CAA) was found 


from 750 to 450 mb, its magnitude decreasing as sur- 
face development increased. 

The composite profiles of LHR are shown in Fig. 
10. The vertical distribution of LHR appears as 
specified in section 3c, with composite maxima at 
700 mb for all development classes. The profiles 
show that the LHR values increased steadily from 
the WK to the SD category. This corresponds to the 
increasing positive near-surface vorticity tenden- 
cies contributed by lath seen in Fig. 7. The small 
negative value contributed by lath seen in Fig. 7 for 
the WK category is a feature that often occurs in the 
finite-difference estimate of the laplacian when the 
cyclone center is located near the boundary of the 
latent heating region, particularly when the latent 
heating is small as is the case here. Although a con- 
sequence of the finite differencing, this result sim- 
ulates the upper-air convergence (and hence surface 
geostrophic vorticity decrease) that would be ex- 
pected in regions adjacent to the upper-air diver- 
gence area attributed to latent heat release. 
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Fig. 11. As in Fig. 8 except for adiabatic heating -cooling (10 * K s' 1 ). 


The composite profiles of dry adiabatic temperature 
change (Sw) are shown in Fig. 1 1. Because of the adi- 
abatic cooling that they represent, these profiles are 
generally of opposite sign to the other thermodynamic 
profiles. In all categories, the magnitude of Su maxi- 
mizes in the upper troposphere at 200 mb. This is a 
result of sufficiently strong upward vertical motions 
(Fig. 12) occurring in the presence of increasing S val- 
ues (Fig. 13). The SD composite shows a secondary 
maximum at 500 mb, resulting from the upward ver- 
tical motion that maximizes at about this level. The 
only positive Soj values occurred in the WK composite 
from 700 to 500 mb, forced by downward vertical mo- 
tions (Fig. 12). The increasing magnitude of the neg- 
ative Suj cooling with increasing surface development 
corresponds to the increasing magnitude of the negati ve 
near-surface geostrophic vorticity tendencies contrib- 
uted by adia seen in Fig. 7. 

Composite vertical motion profiles (Fig. 12) indi- 
cate, as expected, stronger upward vertical motions for 
more intense development categories. The static sta- 


bility profiles (Fig. 13) show hydrostatically stable air 
at all levels with only minor changes with pressure until 
about 300 mb, above which S increases rapidly. Min- 
imum values are seen below 850 mb during the WK, 
WD, and MD categories. However, a small, secondary 
maximum occurs in this layer during strong develop- 
ment. These differences in low-level S values suggest 
that if smaller low-level static stability is a favorable 
condition for cyclone development, then the smaller 
values must exert their influence prior to the onset of 
strong development. Stability values after strong de- 
velopment onset are actually relatively large. This re- 
sult is similar to that reported by Smith and Tsou 
(1988) for a single cyclone case. 

6. Summary and conclusions 

A diagnosis of 12 cyclones that occurred over the 
continental United States and extreme southern Canada 
was performed using near-surface geostrophic vorticity 
tendencies and pressure tendencies derived from an 
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Fig. 12. As in Fig. 8 except for vertical motion (^b s '). 


“extended” form of the Zwack-Okossi (Z-O) equa- 
tion. Pressure tendencies obtained by relaxation of the 
total vorticity tendencies were compared to observed 
sea level pressure changes over the 6-h periods ending 
at each diagnosis time. The comparisons between the 
two fields were very good, as suggested by the average 
correlation coefficient of 0.81 and their comparable 
magnitudes. 

To examine the forcing processes influencing devel- 
opment in the 12 cyclones, composites were generated 
by placing the diagnostic results for each map time into 
one of four categories of development based on the Z- 
O surface pressure tendencies at the grid point of the 
cyclone center. The results show that cyclonic vorticity 
advection was the main contributor to surface devel- 
opment regardless of the cyclone development rate. 
This cyclonic vorticity advection was dominant in the 
upper troposphere and became more prominent as the 
strength of surface development increased. The main 
opposition to development in the weakening and 
weakly developing categories was provided by tem- 


perature advection, as cold-air advection in the lower 
half of the troposphere dominated upper-level warm- 
air advection. Temperature advection contributed to de- 
velopment only during strong development, when 
warm-air advection that maximized in the upper tro- 
posphere was accompanied by a secondary maximum 
in the lower troposphere and the midtropospheric cold- 
air advection decreased in magnitude. Adiabatic pro- 
cesses provided the main opposition to development in 
the periods of moderate and strong cyclone deepening, 
while latent heat release provided increasingly positive 
contributions as the degree of surface development in- 
creased. Only during strong development did vorticity 
advection, temperature advection, and latent heating all 
contribute to development. 

Based on these results we are led to the expected 
conclusion that cyclonic vorticity advection and warm- 
air advection provide upper-level support for cyclogen- 
esis during all phases of cyclone evolution. However, 
somewhat surprisingly, despite the persistent upper- 
level warm-air advection, vertically integrated temper- 
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ature advection supports development only during 
strong development periods. This stands in contrast to 
vertically integrated vorticity advection and latent heat 
release, which are persistent cyclogenetic mechanisms. 
Finally, transitions from weaker to stronger develop- 
ment categories require that the total cyclogenetic 
mechanisms overcome the cyclolytic influence of dry- 
adiabatic cooling. This cooling increases as the cyclone 
intensifies in association with increased upward vertical 
motions and, in later stages, increased static stability. 
The ability to overcome this increased cooling is real- 
ized as 1 ) integrated cyclonic vorticity advection and 
latent heat release increase and 2) upper-level warm- 
air advection is joined by a secondary warm-air advec- 
tion maximum in the lower troposphere, thus negating 
the inhibiting influence of the low to midlevel cold-air 
advection that dominates this term during weaker de- 
velopment stages. 
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